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The problem of resonant tunneling through a quantum dot weakly coupled to spinless Tomonaga- 
Luttinger liquids has been studied. We compute the linear conductance due to sequential tunneling 
processes upon employing a master equation approach. Besides the previously used lowest-order 
golden rule rates describing uncorrelated sequential tunneling (UST) processes, we systematically 
include higher-order correlated sequential tunneling (CST) diagrams within the standard Weisskopf- 
Wigner approximation. We provide estimates for the parameter regions where CST effects can be 
important. Focusing mainly on the temperature dependence of the peak conductance, we discuss 
the relation of these findings to previous theoretical and experimental results. 
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I. INTRODUCTION 

The startling properties of one-dimensional (ID) in- 
teracting electrons, commonly referred to as (non- 
chiral) Tomonaga-Luttinger liquid (TLL) behavior see 
Refs.il 

for reviews, have recently moved into the focus 
of attention of the mesoscopic physics community. This 
was in particular prompted by the successful demonstra- 
tion of electrical transport experiments for a variety of ID 
materials, such as semiconductor quantum wires^ frac- 
tional quantum Hall edge states^ and single-wall carbon 
nanotubes (SWNTs) i In this paper, we mainly focus on 
electrical transport in SWNTs to keep the discussion con- 
crete. Nevertheless, our results apply also to other sys- 
tems. Shortly after the theoretical proposal of TLL be- 
havior in individual metallic SWNTs, 8,9 the first exper- 
imental evidence for this peculiar many-body state was 
reportedii2iAiii2ii£ The expected TLL power-law scaling 
in the energy-dependent tunneling density of states^ 
in SWNTs has been verified experimentally. In later 
experiments^*^ transport through an intrinsic quantum 
dot formed by a double barrier within the SWNT has 
been probed, allowing one to study the well-known reso- 
nant or sequential tunneling including Coulomb blockade 
phenomenafiii2iiS but now for the case of leads formed 
by strongly correlated electrons. When varying an exter- 
nally applied gate voltage, the linear conductance then 
typically displays a sequence of peaks, which can be in- 
terpreted as Coulomb oscillations or resonant tunneling 
peaks, depending on the parameter regime. The peak 
spacing is governed by the charging energy E c and the 
plasmon level spacing e on the dot (the latter coincides 
with the single-particle level spacing E s for the case of 
noninteracting electrons). We mention in passing that 
a similar double-barrier experiment has been performed 
for semiconductor quantum wires, where in fact a TLL 
power-law temperature dependence of the peak conduc- 
tance was reported^ In the nanotube experiments! 1 ^ 6 - 
two intramolecular barriers have been created within an 
individual metallic SWNT, and a power-law tempera- 



ture dependence of the conductance maximum has been 
reported^ However, the exponent turned out to be in- 
consistent with expectations based on the sequential tun- 
neling theory for a TLL; see below. The experiments 
clearly indicate that a detailed theoretical understand- 
ing of transport through a double barrier in a TLL is 
required. 

Since the initial theoretical work on this topic, 14,20 
the double-barrier problem in a TLL has at- 
tracted a significant amount of attention among 
theoristsi 21 i 22 . 23 i 24 i 25 i 26 . 27 i 28 i 29 Those works show that 
the electron-electron interactions present in ID leads 
imply significant deviations from the conventional the- 
ory for Fermi liquid leads^i* 8 - This is most easily seen 
through the change of the temperature dependence of 
the linear conductance peak height, which now exhibits a 
typical TLL power law G max (T) oc T a , where a depends 
on the interaction strength. As is discussed in detail 
below, for a TLL, the exponent a reveals the particular 
tunneling mechanism ruling transport through the dot. 
SWNTs possess an additional flavor degeneracy, 89 
which (together with the spin degree of freedom) can 
be captured by an effective four-channel TLL model. 
For an analysis of the linear conductance on resonance, 
which is the focus of our work, however, a spinless 
single-channel TLL turns out to be sufficient i 1 ^* 2 ^* 2 ^ 
Within the standard bosonization method) 3 " the TLL is 
then characterized by a boson field d{x) with conjugate 
momentum Tl(x). Including a symmetric double barrier, 
composed of short-range scattering centers at ±xo/2 of 
strength Ui mp , the basic Hamiltonian isi- 4 *^ 



Hs 



imp 



(i) 



where vf is the Fermi velocity, 

ffimp = U iiap ^2cos[V4n'&(±xo/2)}, 
± 

and we put h = hs = 1 throughout this paper. Asym- 
metric barriers can also be studied using our approach 
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below, but to keep the discussion as simple as possible, 
we restrict ourselves to the symmetric case alone. The 
electron-electron interaction strength in the leads is mea- 
sured in terms of the standard dimensionless TLL pa- 
rameter g, where g = 1 for a Fermi gas and g < 1 for 
repulsive interactions^* We then have an intrinsic quan- 
tum dot within the TLL, with a plasmon level spacing 
e = E s /g, where E s = nvp/xo is the single-particle level 
spacing, and a charging energy E c = E s /g 2 . (Note that 
our convention for E c differs from the standard onei*^) 
Finally, the coupling to an applied bias voltage V and to 
an external gate voltage Vg, acting onto the dot's elec- 
trons via a capacitance Gg, is encoded in 14,20 



-e(VN/2 + cV G ri), 



(2) 



where c = Gg/(G + Cg) with the island capacitance G. 
In Eq. J2J, we used the notation 



N 



"0, 



(3) 



where — eN gives the charge difference between the left 
and right leads, —en is the total electronic charge occu- 
pying the dot, and no describes a possible offset charge. 
Note that with our conventions, N decreases (increases) 
when electrons are transferred towards the right (left) 
lead. Likewise, n increases (decreases) for tunneling onto 
(out) of the dot. The Hamiltonian Q has been studied 
in most of the previous works on the subject and also 
forms the basis of our work. 

The linear conductance through the dot is periodic in 
-/Vo = ceVG/E s , and we can restrict ourselves to one pe- 
riod, < Nq < 1, where G has a resonance peak at 
A^o = 1/2. One can compute G analytically for the non- 
interacting limit, g = 1, by refermionizing this model, 28 



G 



Go 



dE 



w 



4Tcosh 2 (E/2T) cos 2 {ir[N + E/E s ]) 



(4) 

with G = e 2 /h and w = (4 - A 2 ) 2 /[8A(4 + A 2 )], where 
A = 7r£/i mp /D for the bandwidth D, which provides a 
high-energy cut-off. For strong barriers, Eq. leads to 
the standard Breit-Wigner resonant tunneling line shape 
with linewidth Tq = wE s /n. (For g = 1, the infinite- 
barrier limit is reached already for A = 2, where the 
associated phase shift is in the unitary limit.) Unfortu- 
nately, as the model Q is not integrable for g < 1, ex- 
act solutions covering a wide parameter range of interest 
for this transport problem are out of reach. Analytical 
progress then generally has to rely on approximations. 
Initial workia^S*Si pursued perturbative approaches, us- 
ing the renormalization group, instanton methods or cu- 
mulant expansions, both for strong and weak barriers. 
Furusaki 2 * presented a detailed study of the uncorrelated 
sequential tunneling (UST) regime, including also cotun- 
neling contributions important away from the resonance 



peak. This regime allows for a master equation approach, 
whose validity requires that the barriers are sufficiently 
strong and that the temperature T is sufficiently high. 
For T well above e but still below E c , the discrete nature 
of the plasmon modes on the dot is not relevant while 
Coulomb blockade still exists. Then the corresponding 
results of Ref . ITU apply. In particular, one obtains power- 
law scaling for the peak conductance with a = 2(1/ g— 1), 
implying that in this case the double barrier effectively 
acts as a single impurity. In what follows, we only discuss 
the case T < e. 

Keeping only rates to lowest order in Tq from the TLL 
leads onto the island (dot), which is equivalent to taking 
the standard UST mechanism, the conductance is given 
byS2. 



Gust T {wT/D) 



-l+l/g 



Go 



4r(l/ 5 )Tcosh(<5/2T) 



1 

2~9 



iS 
2nT 



(5) 



where T(x) denotes the Gamma function, and 5 = 
E c \Nq — 1/2 1 measures the distance from the resonance 
peak. The hybridization Tq = 2ttpq/S, 2 can be expressed 
in terms of the ID density of states, p$ = 1/(ttvf), and 
the hopping amplitude A onto the dotjSi 

A = Tr-^l + 1/5)^(1 +g)(nU imp /D)- 1 /aD. (6) 

The line shape 10 is rather close to the Fermi-liquid form 
(|4"|) for large barrier heights and evidently characterized 
by a linear T dependence of the linewidth. In the tails of 
a peak, the conductance Q vanishes exponentially, but 
then also (elastic) cotunneling has to be included ^LSLSSi 
Note that the UST peak conductance in Eq. 10 scales 
as G max cx T~ 2+1 / g , so that the power-law exponent is 
a uST = — 2 + l/g. Finally, at low temperatures, instead 
of sequential tunneling, coherent resonant tunneling is 
possible, characterized by non-Lorentzian universal line 
shapes,— where the linewidth scales as T 1 ^ 9 . In this pa- 
per, we only address the incoherent regime by employing 
a Markovian master equation approach. Moreover, we 
focus on the temperature dependence of the peak con- 
ductance, for which cotunneling is always a subleading 
process^ Therefore we neglect cotunneling throughout 
this work. This is not a fundamental restriction to our 
approach, but implies some technical simplifications. On 
resonance, the basically only restriction for the validity 
of our master equation approach is then given by the 
condition G max <GoS 

Remarkably, the available experimental data obtained 
in SWNTsi* seem incompatible with Eq. JSJ, since the 
observed temperature dependence of the conductance 
peak height G max does not follow the UST scaling but 
rather suggests a T~ 3+2 / 9 power law. The exponent 
&cst = S + 2/g has been proposed to arise from a 
"correlated sequential tunneling" (CST) mechanism. 24 
CST processes are sequential tunneling processes that 
cannot be subdivided into two uncorrelated steps, yet 
still can be captured in a master equation framework. 
Such processes have been studied previously in the chem- 
ical physics community, e.g., in the context of electron 
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transfer reactions through an intermediate bridge state. 
In particular, Hu and Mukamel— have treated the se- 
quential tunneling regime for this problem using a very 
similar master equation theory. Here we compute the se- 
quential tunneling current through a double barrier in a 
TLL, see Eq. JJJ, beyond the lowest order in Tq, imply- 
ing modifications to Eq. We analyze to what extent 
such processes could indeed cause power-law scaling in 
G max (T) with the CST exponent cucst- 

The remainder of the paper is structured as follows. 
In Sec. [H]we introduce the general master equation ap- 
proach, and apply it to the regime of linear transport. 
The transition rates entering the master equation are 
evaluated in Sec. IIIII and in Sec. IIVI explicit results for 
the temperature dependence of the conductance peak, 
G maK (T), are presented. We conclude by discussing our 
results and their relation to other (experimental as well 
as theoretical) work in Sec.[V] Technical details concern- 
ing Sec. IIIII can been found in the Appendix. 



II. MARKOVIAN MASTER EQUATION 

Let us start with the case of large tunneling barriers 
EAmp, which can be describe d within a dual version of 
Eq. CP, seeRefs, IT320l 22 30. In this regime, the dynam- 
ics is dominated by tunneling events connecting minima 
of the periodic potential H lmv = 2£/i mp cos(7rTV) cos(7rn) 
in the (TV, n)-plane^ Such tunneling events induce a 
change n — > n± 1 for tunneling onto/out of the dot, and 
TV — » TV =F 1 for tunneling towards the right/left. Hence 
transfer of one unit of charge across the complete dou- 
ble barrier structure requires TV — ► TV ± 2. The current 
through the double barrier is then given by 



I=- 2 (N), 



(7) 



where the expectation value stands for a quantum- 
statistical average with Hamiltonian , and Eq. J7J has 
to be evaluated in the stationary long-time limit. The 
discrete dynamics underlying Eq. J7J can be described 
by a master equation^ for the probability Pjv(n, i) of 
being in state (TV, n) at time t. 

Master equations have previously been employed for 
the non interacting casejiiiSiiS and for the TLL case 
to lowest order in Tom^^ The rates entering the mas- 
ter equation can be extracted as irreducible diagrams 
for the self-energyjSi which to lowest order are simple 
golden rule rates. These first-order UST contributions 
to the transition rates imply TV — > TV ± 1 jumps. Below 
we will systematically take into account transition rates 
up to second order in Tq. It turns out that the second- 
order contributions to those rates are plagued by nontriv- 
ial divergences, which require a resummation of higher- 
order diagrams. This resummation is done below by em- 
ploying the Weisskopf-Wigner approximation^ On the 
one hand, this procedure yields a direct CST process, 
TV — > TV ± 2. On the other hand, additional (indirect) 



CST contributions to the transition rates for TV — > TV± 1 
arise. Apart from the need of a regularization scheme, 
the same situation is encountered in the theory of bridged 
electron transfer reactions^ 

The master equation for P/v(n, t) then has the general 
structure^ 

P N (n,t) = -j{n)P N (n,t) (8) 
+T f R (n + l)P N+1 (n + l,t) + T b L {n + l)Pv-i(n + l,t) 
+r{(n - l)P N+l {n - l,t) + T b R {n - l)Pv-i(n - l,t) 
+T f CST (n)P N+2 (n,t) + T b CST (n)P N ^ 2 {n,t) 
-T f CST (n)P N (n, t) - r b CST (n)P N (n, t), 

where r£y^(n) denotes the forward/backward rate for a 
transition over the left/right barrier, having started with 

n electrons on the dot. In addition, T^Jg T {n) denotes 
the forward/backward rate for a direct CST transition 
TV — > N ± 2. Moreover, we use 



7 (n) = r£(n) + r{(n) + T b R (n) + T» L (n), 



(9) 



which is related to the linewidth of the state (TV, n) . It 
is useful to also introduce the probability 



o(n,t)= J2 p N(n,t) 



(10) 



N=-c 



for the dot being occupied with n electrons at time t. In 
order to calculate the current /, we insert Eq. JSJ into 
Eq. J7J, which yields after some algebra 



i = -- y 

2 ^ 

n — — oo 



[T R (n)+T f L (n)-T b R {n)-T b L (n) 



2L 



/ 

CST 



(n) 



2T b CST (n) p(n,t-> oo). 



(11) 



Summing both sides of Eq. JSJ over TV, we obtain a mas- 
ter equation for p(n, t) directly, 



p(n,t) = —"/(n)p(n,t) 

Il(n+l)+r£(n. 

r{(n-l)+r^(n- 



(12) 



p(n + l,t) 
p(n ~ 1, t). 



The stationary solution p{n) — p(n,t — > oo) follows by 
requiring p(n) = 0, which yields the detailed balance 
relation 



p(n) _ r^,(n + l)+r^(ra + l) 
p(n + l)~ r f L (n)+r b R {n) 



(13) 



Taking into account conservation of probability, 
S n =-oo P( n ) ~ 1' tbis relation can be solved recursively. 
We note that the direct CST rates r£/^ T (n) do not 
appear in Eq. (|12f> . since they do not alter the net 
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population of the dot. However, they do appear in the 
current 

Let us now focus on the linear transport regime and 
sufficiently low temperatures, eV, T <C E Cl where at most 
two charge states n are allowed on the dot due to charging 
effects. 19 Without loss of generality, we may choose n = 
and n = — 1 to label those states. In the linear transport 
regime, we can disregard the rates r{(0), 1), r^(0), 
and T b L (—l), since they involve energetically forbidden 
states with dot occupation numbers n = +1 or n = —2. 
The recursive solution of the detailed balance relation 
(|T3j) then yields 



Combining Eqs. i|14|) and I jllj l. we obtain / 
with the standard contributio n 1 !) 1 ?! 1 ?! 22 ! 2 ? 



h + h 



h 



r£(o)r£(-i)-rS 



-i)rl(o) 



r£(o) + r£(-i) + i*(-i) + r£(o)' 



(15) 



p(0) 

K-i) 



r((-i) 



i) 



r£(o) + r£(o) + r£(-i) 



(14) 



r fc L (o) + r^(o) 



r h L (o) + r^(o) + ri(-i) + r^(-i) 



and an additional term caused by direct CST rates, 



CST 



(0) 



r{(-i) + r^(-i) 



+ 



CST 



(-1) 



rksr(-i) r^(o) + 4(o) 



ri(o) + r£(-i) + i*(-i) 



r h L (0) 



(16) 



So far we have discussed a general procedure to deter- second order in To, 
mine the current in the linear regime. To make progress, 

the transition rates entering the above equations must be T^(n) = L^ 1 ^ (n) + T^ 2 ' (n), (17) 

computed for the Hamiltonian in Eq. (Q. 

where we keep only the n = 0,-1 states. 



III. TRANSITION RATES 

We now systematically compute all rates entering 
Eqs. (|T5|) and (|TT)|l . Rates of order higher than T 2 , will 
be included approximately within the Weisskopf-Wigncr 
scheme. For consistency, besides the direct CST rates, it 
is also mandatory to include all indirect CST contribu- 
tions to the rates T^(n) with A = L/R and v = f /b. All 
these transition rates can be extracted as the irreducible 
components of an exact perturbation series expression 
for the probability distribution P/v(n,i). The latter cor- 
responds to the diagonal element of the reduced density 
matrix (RDM), which in turn allows for a path-integral 
representation^ Tracing out the Gaussian TLL modes 
away from the barriers at x = ±xo/2, one obtains an ef- 
fective action which is equivalent to the action of a quan- 
tum Brownian particle hopping in the (N, n)-plane^Si 
Then a path can be visualized in the (N, 7V')-plane of 
the RDM, see Fig. ^ with a corresponding dynamics 
in the (n, n')-plane. Alternatively, the rate expressions 
given below can also be derived using the projection op- 
erator formalismi£ii22i2i The rates T^(n) for a transition 
N — > N ± 1 then have contributions of first and at least 



A. First-order rates 

The first-order contribution is well-known) 14 i 20 i 21 i 22 i 30 
and schematically depicted in Fig. [^a). There are two 
independent (uncorrelated) steps, one from the left TLL 
"lead" onto the island, and another to the right TLL 
"lead". Using the hopping matrix clement A in Eq. J()J), 
these steps individually correspond to irreducible golden 
rule rates, 



r^' (1) (n) = ^Re 



A 2 



2 jo 



where W^(t) 
functions 21 



dtwcp[iE Xv {n)t-Wv{t)}, (18) 
W + (t) + W_(i), with the correlation 

J±(w) 



/•DC 

W±{t) = du 
Jo 



-{[1 - cos(wi)] (19) 



x coth(w/2T) -I- i sin(wt)}. 
The spectral densities follow as 

j/D 



cue 



2.9 



l + 2eY^ 5(w - M±(m)e) 



(20) 
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where M + (m) = 2m— 1, M—(m) — 2m, and e = i? s /<7 = 
irvp/(gxo) is the plasmon level spacing on the dot. The 
(5-peaks are a result of the finite level spacing on the dot, 
while the first part in Eq. (|2(J|) yields a standard Ohmic 
spectral density. The bandwidth D is taken as smooth 
(exponential) ultraviolet cutoff for the model 0). These 
correlation functions arise in the process of integrating 
out the TLL modes away from the barrier, see Ref. 30 for 
a detailed review of this procedure. Finally, the energies 
E\ tV {n) appearing in Eq. ljT%|) are defined as 

E Rf {n + l) = -E Rb (n) = v{n + l)-eV/2, 
E Lb (n+l) = -E Lf (n)=ti(n + l) + eV/2, (21) 

with the electrochemical potential 

fjb(n + 1) = 2E c (n - n - ecV G /2E c + 1/2). 

The correlation function Ws(t) = W+ + W- = Wohm + 
Wdot can be decomposed into two different parts ; 21 i 2 ? 
namely an Ohmic part Wohm(t) and an oscillatory part 
Wdot(i). Here the first contribution comes from the 
smooth part in J±(u>), and takes the standard form 

Wohm(t) = ~ln[(D/T)smh\%Tt\]+i(Tr/g)sgn(t) 
9 

= S h m (t)+iRohm(t). (22) 

At low temperatures, T<e, the dot correlation function 
is given by its T — limiting form 

^dot(*) = ~ln( — j, (23) 

up to exponentially small corrections in y = e~ E / T . This 
function arises due to the finite level spacing, and is pe- 
riodic in t with period r e = 2ir/e. It can therefore be 
expanded in a Fourier series? 2 ^ leading to 

F\ (1) (n)= d p {e)T 0iim {Ex„(n)-pe). (24) 

p— — oo 

The Fourier coefficients d p (e) can be found in the Ap- 
pendix, and from Eq. Q22JI. one gets^l 

A 2 e E / 2T ( D V- 1/9 (I iE\ 2 

To ^ E '-WrOJi){2^f) T \Yg + 2^T) ' 

(25) 

which in turn directly leads to Eq. We note that the 
first-order forward/backward rates (|18|l fulfill a detailed 
balance relation 

(n - 1) = e -^/(»)/^r{' (1) (n), (26) 

which formally follows from the reflection property 
W±(-t) = Wl(t) = W±(t- i/T) of the above corre- 
lation functions. 
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TABLE I: All transition processes of order Fq contributing to 
the master equation. The off-diagonal states a, b and c are 
specified in Fig. 



B. Indirect CST rates 

The Tq contributions to the rate for TV — > TV — 1, 
which we shall call "indirect transitions" , require a care- 
ful counting of all possible transitions consisting of four 
jumps in the (TV, TV')-plane, see Figure and Tabled 
There are eight irreducible diagrams for the forward rate, 
plus their complex conjugates, which can however be in- 
cluded by taking twice the real part. Similarly, there 
are eight diagrams for the backward rate. Each irre- 
ducible second-order diagram then gives a triple integral 
over three times Ti,T2 and T3, which represent the time 
spent in the corresponding state of the reduced density 
matrix (RDM)! 2 ^*^ To give a concrete example, the dia- 
gram denoted by k = 7 in Table [I] is drawn schematically 
in Fig. |5Ja). In particular, since we neglect cotunneling 
processes, after two jumps we are always in a diagonal 
state of the RDM, amounting to a real (as opposed to vir- 
tual) occupation of the corresponding state. Therefore T2 
in those expressions will always have the meaning of the 
time spent in the respective intermediate diagonal state. 
Similarly, for the "direct" CST rate with TV ->■ TV — 2, 
there are two such triple-integral contributions (plus their 
complex conjugate diagrams). At this point, we note that 
we label the direct transitions as CST1/2 and the indi- 
rect transitions with the index k — 1, ... ,8, although in 
principle all second-order processes taken into account 
here are correlated sequential tunneling processes. 

An important property of the irreducible second-order 
diagrams is that for finite plasmon level spacing e, each 
of them contains a divergence. This is not a trivial diver- 
gence as we work with irreducible diagrams. Formally, 
this infrared divergence comes from the T2 integration 
that extends from to 00 but has an integrand which 
is ultimately periodic in t%. This implies that one has 
to effectively include higher-order diagrams, which is of 
course impossible to achieve in an exact way. Here we 
use the Weisskopf-Wigner approximation^ to regularize 
the second-order diagrams. Physically, the intermediate 
state has a finite lifetime linked to the linewidth param- 
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eter 7(71) in Eq. JJjJ. We then introduce a factor e _7( -™) T2 
for the r 2 -integrations, where the linewidth 7(71) has to 
be computed self-consistently by requiring that Eq. 
holds. This linewidth 7(71) then enters the master equa- 
tions JHJl and |Q . 

The eight regularized rates r^'^ are given explicitly 
in Table [n] At low temperatures, again the T = form 



1 /l + e -( £ /- D + i£ *) 
W A (t) = W+(t) - W-(t) = --In 1 + e - E/g 



(27) 



holds up to exponentially small corrections of order y 
e~ e l T , and 



8 

fe=i 



",(2) 



(28) 



This rate depends on 7(71), which needs to be determined 
self-consistently. Note that W&(t) is also periodic in t 
with period r e . 

From now on, we always focus on the resonance peak, 
where all energies can be put to E\ v {n) = 0. In ad- 
dition, since on resonance /i(— 1) = /i(0), we also have 

7(— 1) = 7(0) = 7. The calculation of is com- 

paratively simple for the diagrams k ~ 1,2,7,8, since 
only the periodic correlation function appears in the 
respective bracketed terms, which allows for a Fourier se- 
ries expansion. Then the triple integrals factorize, and 
the quickly converging Fourier sums can easily be per- 
formed numerically. For the benefit of the interested 
reader, in order to illustrate this procedure, the evalu- 
ation of diagrams of this first class is discussed in full 
detail in the Appendix. The remaining diagrams of the 
second class (k — 3, 4, 5, 6) are more difficult to han- 
dle, since the bracketed terms now involve the correla- 
tion function Wj. We have therefore evaluated the re- 
spective triple integrals numerically as a function of the 
linewidth 7. This can be done either via direct numeri- 
cal quadrature (trapezoidal rule), or using Monte Carlo 
integration^. The latter approach is more suitable for 
large 7, where stochastic error bars can be made ar- 
bitrarily small with only modest computational effort. 
Fortunately, the 7-dependence of the four diagrams of 
the second class (k = 3,4,5,6) turns out to be identical 
to the one of the first class, which is given by Eq. I|29|) 
below, see also Eq. (|A3ll in the Appendix. A fit of the 
numerical results for several 7 then allows to accurately 
extract the parameters A g ^,B g ^ and C g ^ for those di- 
agrams as well. Finally, we summarize all indirect CST 
contributions as 



W2) 



A 4 



-A n 



7 



B + C, 



(29) 



where dimensionless parameters A g ,B g , and C g follow 
by summing over the respective values A g k, B g k,C g ^ 
for these eight diagrams (including their complex conju- 
gates), see AppendixlAl These parameters depend on the 



TLL parameter g and on the dimensionless temperature 
T/s. 

Of primary interest is then the temperature depen- 
dence of the self-consistently determined linewidth pa- 
rameter 7. Unfortunately, as we discussed above, it seems 
impossible to evaluate A g ,B g , and C g analytically, even 
in the non- interacting limit g = 1. However, numerically 
we can obtain them for given (g,T), see Appendix lAl for 

Since the master 



details, and we find B„ 



Ag < C g . 



equation approach holds only for 7 <C e, it is clear that 
the B g term in Eq. I|29[) can be neglected for all practical 
purposes. Numerical results for A g and C g for different 
g and T are shown in Table Hill and Fig. 21 respectively. 
Equation (|29|l indicates that C g follows from the large- 
7 behavior of T v ^ 2 \ while determining A g requires the 
small-7 solution of the relevant triple integral. Unfor- 
tunately, the latter is numerically rather expensive at 
low temperatures, and hence we can specify A g only for 
moderately low T, see Table Hill while Fig. 0] covers our 
results for C g down to T = O.Ole. Evidently, the tem- 
perature dependence of C g becomes rather weak at low 
temperatures, C g (T) ~ const., an observation supported 
by analytical arguments given in Appendix [S] Based on 
these arguments, we expect that even for g = 0.9, where 
Fig. ^suggests a significant T-dependence, at sufficiently 
low T the quantity C g becomes independent of tempera- 
ture. 

We mention in passing that in contrast to the first- 
order contributions in Sec. 1111 Al the above indirect CST 
contributions do not obey detailed balance. This is true 
although the correlation functions entering these rates in 
Table ITTl still have the reflection property. The violation 
of the detailed balance relation follows directly by inspec- 
tion of the triple integrals in Table ITTl i- e -i forward and 
backward rates are not linked by a relation of the form 
(|26|l . Of course, the total rates and populations of the 
states are still linked by detailed balance, but there is 
no reason why individual rates should obey Eq. 1)26(1. A 
simple example for this fact is given by superexchange 
rates in electron transfer theory^2i2i see also Ref. for 
related observations. 



C. Direct CST rates 

The direct CST rates T^j ST (n) for a transition N — > 
N — 2 have only contributions of at least order Tq. There 
are two possible transition processes CST1 and CST2, 
see Table UJ and Figs. ^ an d EJb). The corresponding 
transition rates are simply 



r 



f ( 

CST,1\ 



,/,(2) 
R,7 



(n), 



CST.2 



(n) 



R,8 



(n). 



These relations hold not only on resonance but in gen- 
eral. Since CST transitions do not alter the stationary 
population p(n), they do not enter the linewidth The 
rates r^ ST (n) per se are also not subject to a detailed 
balance relation l|26|l . Note that the cotunneling diagram 
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Forward N —* N — 1 rates of (at least) order F 2 , over the right barrier: 
r^(n) = -2^Re jf dr, 

x e i[-Sz,bCn)ri+BH / (n)T3] e -W^(T 1 )-WE(T 3 ) r e Wi(T 1 +T 2 )-W^(r 2 )-W^(r 1 +r 2 + r 3 ) + W^(T 2 +T3) _ jl g- 7 (n)T 2 



r^(n) = -2^Re ,f »dr, 
r^ 2) H = -2^Re f °°dn 

X e »[-BH/(n)Ti+B H/ (n)r 3 ] e -W s (T 1 )-W I; (r3) Lw E (-n+-r 2 )- W s (r 2 ) - W E (ri+T 2 +T 3 )+Ws fa +t 3 ) _ jl g-7(n)T 2 

4' ( 4 2) H = -2#Re /-dTi 

X e »[-BH/(n)Ti-B H/ (n)T 3 ] e -W s (T 1 )-W*(T3) F- W s fa +r 2 ) + (r 2 ) + W s fa +r 2 + T 3 ) - W s (r 2 + T 3 ) _ ]] g- 7 (n)T 2 

r^(n) = -2^Re J^dr, 

x e i[ B B/(«)'-l- B B6(™-l)T3] e -'»'ECri)-W^g(r3) jgWs fa+r 2 ) - W£ fa )- W s fa +r 2 +r 3 ) + W£ (t 2 +t 3 ) j_j g — y(n)T 2 

r« 2) (n) = -2^Re / °° dr. 

x e l [ B H/(™M+EH6(n-l)T 3 ] e -W E ( T1 )-W s (T 3 ) ^ e -M , E (T 1 +T 2 ) + W£(r 2 ) + VK s (r 1 +T 2 +T 3 )-W£(r 2 +T 3 ) _ j] g-TOfa 

r^ 7 2) (n) = -2^Re J^dr, 

x e i[ B B/(«)-ri+- s; £/(«- 1 )T3] e -W r s (ri)-M' E (-r3) jgW^ fa+T 2 )-W A (t 2 )- W A (n +t 2 +t 3 ) + W a fa+r 3 ) _ j] e -T(«)-r 2 

r^ 2) H = -2^Re J^dr, 

x e l [- B H/(™)- r i-- B l,/(' l - 1 ) T 3] e -VK s (T 1 )-W*(T 3 ) ^ e -WA(Ti+T 2 )+W A (r 2 )+WA(Ti+T 2 +T3)-W A (T 2 +T 3 ) _ e -TC«)T2 



TABLE II: The 8 irreducible rate expressions of at least order Fq, corresponding to the forward-rate diagrams through the 
right barrier, F^;' 2 ', see Eq. I|280 and Table U] The corresponding rates through the left barrier, T^f\n — 1), follow by 
substituting E R f/ Lb (n) — > E L f/ Rb (n — 1) and E L f/ Rb (n — 1) — » E R j / Lb (n). The backward rates V b R ? (n — 1) and r^'^^n) 

can be obtained from the forward rates r^ 2 '(n) and r£'^ 2 '(n — 1) using the substitutions E Lb / R j(n) — > E L j/ Rh (n — 1) and 
E Lf/Rb {n - 1) -> E Lb/Rf (n), respectively. 





T = 0.1e 


T = 0.2e 


T = 0.5e 


5 = 0.4 


3.37 x 1(T 8 


-1.56 x 10~ B 


-1.12 x 10~ B 


5 = 0.5 


3.00 x 10" 8 


-2.02 x 10~ B 


-6.07 x 10 _s 


5 = 0.6 


2.35 x 10~ s 


9.71 x 10~ a 


2.83 x 10 _Y 


5 = 0.7 


-1.97 x 10~ 8 


2.32 x 10~ s 


1.34 x 10~ b 


5 = 0.8 


-1.42 x 10 w 


7.68 x 10~ 8 


3.86 x 10~ b 


5 = 0.9 


-3.76 x 10 _Y 


7.35 x 10"' 


7.85 x 10~ b 



TABLE III: Numerical results for the dimensionless parameter A g in Eq. 129H for various g and T at D — We. 



in Fig. IHb) is subleading on resonance and therefore not Neglecting the _B s -term as explained above, the self- 
taken into account here. consistency equation for 7 can then be written as 



D. Linewidth 

On resonance, E\ u (n) = 0, and hence forward and 
backward rates are equal, 



r /,(l)/(2) _ p6,(l)/(2) _ p(l)/(2) 



(31) 



7 = 4(r« + r( 2 )) = 4r« + ^ (-~A g + , 

£ J \ 7 



C g . (32) 



This quadratic equation has two real solutions, as long 
as the dimensionless parameter £ < 1, where 



( e 3 r (i) + A 4 c ff ) 2 ' 



(33) 
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If indeed £ < 1, these solutions are given by 

7± = 4 



r (i) , A4g 3 



For ^ < 1, we thus have 
7 = 7+ 



(34) 



(35) 



and 7_ = £7+/4. For £ > 1, the two solutions are com- 
plex valued, and 7 acquires an imaginary part. If this 
happens, the description in terms of the master equation 
in combination with the Weisskopf-Wigner approxima- 
tion is questionable and not used below. We identify the 
linewidth with the solution 7+ , since we recover the well- 
known result 20 7 = 4rW when second-order rates are 
neglected, while then 7_ = 0. The requirement £ < 1 
typically results in a temperature 7} determined by 



r«(T}) + ^ 



below which our Weisskopf-Wigner theory becomes prob- 
lematic. This equation yields Ti provided A g (T) and C g 
(assumed independent of temperature) are known. How- 
ever, since no reliable low-T estimates for A g are avail- 
able, see Table IIIII it is often difficult to provide good 
estimates for 7]. Of course, the validity of the master 
equation in addition always requires G max <C Go- In 
what follows, the parameter C g is assumed to be temper- 
ature independent and given by its value at T/e = 0.01 
in Fig.m Although we expect C g {T) to be constant (see 
above), we cannot rule out that this is approximative. 

The linewidth j(T) now consists of two contributions. 
The first term in Eq. I|35|) is cx T 1 / 9-1 , while the second 
term is constant due to the T independence of C g . This 
implies a crossover from power-law scaling of 7(T) to a 
basically constant 7 as the temperature is lowered. This 
crossover depends in an essential way on the tunneling 
amplitude A. Results for j(T) from Eq. Ij35[) are shown 
in Fig. Eta), taking A = 6e and D = lOe. Figure Efb) 
shows (for several T) that the validity condition £ < 1, 
with £ defined in Eq. I|33|l . is indeed fulfilled. This choice 
for A reflects rather transparent barriers, where CST ef- 
fects are clearly pronounced. Obviously, at low T, Fig. [3] 
suggests that 7 is essentially independent of temperature. 
This behavior is most pronounced for strong interactions 
(small g). Going towards less transparent barriers, j(T) 
is shown in Fig. for A = 3e. In Sec. [V] we argue 
that the hoppingamplitudes A appropriate for the ex- 
periment in Ref. [l5| and for the numerical simulations in 
Ref. are comparable to this value. Now a crossover 
from the UST power law at high temperatures to an ap- 
proximately T-independent behavior at low T becomes 
apparent. For g < 1/2, we now find £ > 1, implying 
that our approach breaks down for such interactions. Fi- 
nally, for very high barriers, A — > 0, the linewidth is 
aiways dominated by the UST term, in accordance with 



standard reasoning, 2 ?* This is illustrated in Fig. Q where 
results for 7 at A = O.le are depicted. In this regime, 
higher-order corrections are obviously negligible. 

Remarkably, for weak interactions, the linewidth 7 is 
always dominated by the UST result over the entire range 
of temperatures where the master equation approach is 
valid (£ < 1). CST effects then apparently do not have 
a finite domain of observability in the limit of weak in- 
teractions. We can therefore rationalize that previous 
calculations that essentially expand around the noninter- 



acting casi 



,25.26,29 



do not observe a clear CST power-law 



scaling. The fate of CST effects near the noninteract- 
ing limit g = 1 will be further discussed at the end of 
Sec. IIVI We conclude that for CST effects to be observ- 
able, it is essential to allow for rather transparent bar- 
riers, a finite level spacing, and intermediate-to-strong 
interactions. The parameter regime where CST plays a 
role is therefore rather narrow. 



IV. TEMPERATURE DEPENDENCE OF 
CONDUCTANCE PEAK 

The linear conductance G follows directly from 
Eqs. H15|) and (|16[) by performing the derivative with 
respect to the transport voltage V. Moreover, since 
we are interested in the conductance maximum at res- 
onance, we can put /i(n) = u.(n + 1) = 0. For the first- 
order rates, the detailed balance relation l|26[l can be ex- 
ploited. For the second-order rates, we find the relation 
dTf^/dV = -dT b 'W/dV. Finally setting V = 0, we 
obtain 



G 

max 

(T)=G 



max, A 



G 



max.-B ? 



(36) 



where 



G n 



T 7 



(*»)■ , 



G max . B = - e r( 2 )'-2 e r' CST + ^r( 1 )r( 2 ). (37) 

1 7 

Here the prime denotes the derivative d/dV, taken at 
V = 0, and we omit the arguments of the rates, since 
on resonance they are all equal. The dominant term is 
the first term in G max ,A- ln order to show that the sec- 
ond term G max ,B is always negligible against G maXjJ 4, it 
is instructive to evaluate the second class of diagrams 
(k = 3,4,5,6) in Table ITU within a simple approxima- 
tion. Since at long times the Ohmic part in Ws (t) cancels 
out in the square-bracketed terms of the corresponding 
rate expressions in Table ITU one may replace Ws(t) by 
Wdot(i) in those square brackets. After this replacement, 
the Fourier expansion method discussed in the Appendix 
applies to all eight diagrams, with Fourier sums includ- 
ing terms like dK R / L (E)/dE. These sums can easily be 
calculated numerically. For all cases that are described 
below (and many more not shown here), we find that 
G max .s is numerically zero. Although it looks like the 
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direct CST rates Tqst have effectively no influence in 
the end, this is not true since they exactly cancel certain 
contributions to the conductance coming from the indi- 
rect rate r( 2 ). In that sense, the direct CST diagram, 
jumping by two steps along the diagonal of the reduced 
density matrix via an intermediate diagonal state, see 
Fig. ^ but without cutting the diagram in the interme- 
diate state, is crucial for CST effects in the conductance 
maximum G max . Such diagrams are known to cause im- 
portant effects in other systems^ but were previously 
not taken into account since the main focus was on the 
limit A -> 0. 

Hence we find for the conductance maximum 



Gmax (-^J 



(38) 



Judging from our numerical results for G maXj B) the '~' 
should in fact be replaced by an exact equality, although 
we have no analytical proof for this statement. For T > 
T c , we find 7 cx T 1 / 9-1 , leading to the UST result ©. 
However, for T < T Cl j(T) stays approximately constant, 
and hence an approximate power-law behavior follows, 



G n 



oc T 



2/3-3 



(39) 



with the CST exponent «cst = — 3 + 2/g. We stress that 
Eq. (|39l) is not meant in the sense of universal power- 
law scaling behavior. The crossover between UST and 
CST-dominated regimes occurs around a temperature T c 
discussed below. The effective doubling in the exponent 
reflects the physics of this correlated transport process. 
For T <C £, plasmon modes excited on the island will 
correlate electrons in both leads. Since each lead has 
a end-tunneling density of states oc E 1 / 9 ^ 1 , correlated 
transport leads to an effective doubling in the exponent 
due to the presence of two leads. 2 - Results for G max (T) 
at A = 3e are shown in Fig. [S]and follow Eq. 1|39[) at low 
T. For g > 0.8, the master equation breaks down (G max 
exceeds Go), while for g < 1/2, the validity condition 
£ < 1 is violated. Nevertheless, there is a well-defined 
region of applicability, where CST effects are important 
and observable. Finally, for A = O.le, the expected UST 
scaling is recovered, see the inset of Fig. |HJ 

The crossover between the UST and CST regimes is 
characterized by a temperature T c — T C (A,<?), which in 
turn follows from the condition that both contributions 
to 7 in Eq. be equal, T^(T C ) = A 4 C g /e 3 . Although 
rW contains the sum over all Fourier modes, see Eq. Q24[l. 
the zero mode p = dominates on resonance, implying 
the condition doFohm(O) = A 4 G g /e 3 . This leads to the 
crossover temperature 



A 2 G g 



9/(1-9) 



(40) 



where we use the abbreviation 

(l-e- £ /- D ) 1 /9r 2 (l/2 3 ) 



ADT{\/g) 



(2tt/Z)) 1 /s- 1 . 



In Fig. |21 results for the crossover temperature T c are 
shown as a function of A/e for the TLL parameters 
g = 0.6 and g = 0.7, always respecting the validity con- 
ditions £(T C ) < 1 and G max (T c ) < G . Generally, T c 
increases when increasing A and/or decreasing g, i.e., 
for more transparent barriers and/or stronger interac- 
tions. Apparently, for weak interactions, g < 1, the A 
dependence of T c becomes very steep, restricting the CST 
regime to extremely low temperatures for reasonable A. 
In practice, at such low temperatures coherent resonant 
tunneling processes dominate, rendering CST effects un- 
observable. For stronger interactions, however, CST ef- 
fects can be pronounced even for moderately transparent 
barriers at low temperatures. 



V. DISCUSSION 

By using a master equation approach, linear trans- 
port in a TLL with two tunneling barriers forming a 
quantum dot has been studied. We find an approximate 
power-law temperature dependence of the peak conduc- 
tance in the linear transport regime, with a characteristic 
^-dependent exponent, where g is the TLL parameter. 
By including second-order contributions to the tunnel- 
ing rates in combination with a self-consistent Weisskopf- 
Wigner regularization, a comprehensive picture has been 
obtained. For temperatures below the plasmon level 
spacing e = nvp / (gxo) for dot size xq, the dominant 
transport mechanism depends on the transparency A of 
the barriers. For sufficiently transparent barriers, A > e, 
we find that correlated sequential tunneling (CST) is im- 
portant, leading to an approximate power law with ex- 
ponent acsT = — 3 + 2/g, while for very high barriers 
(A — > 0), the uncorrelated sequential tunneling (UST) 
exponent «ust = — 2 + 1/g is recovered. Note that de- 
spite the large A necessary to reach the CST regime, the 
master equation should still apply as the peak conduc- 
tance remains small. We have determined the crossover 
temperature T c separating both regimes. As a function 
of the physical parameters (g, T, A, e, D), it is given by 
Eq. (|40l) , where the dimensionless number C g can be read 
off from Fig. 2| For T < T c , the CST mechanism is effec- 
tive, while for T > T c , the UST picture is recovered. The 
peak conductance in the CST regime has been derived 
and is given by Eq. I|38|l , with the well-known first-order 
rates specified in Eq. I|24|l . The linewidth parame- 
ter 7 as a function of the physical parameters is given 
in Eq. 1)35(1 . These relations allow for a comparison to 
existing work and provide estimates for the parameter 
regime where such effects should be of importance. No- 
tice that the modification of the standard picture of se- 
quential tunneling here arises due to a renormalization of 
7 by higher-order processes, and we have given a physical 
explanation for this mechanism above. 

Our findings regarding approximate CST power laws 
in G max (T) are consistent with recent numerically exact 
real-time quantum Monte Carlo simulations 28 and also 
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with experimental observations. Let us first discuss the 
experimental work on SWNTs^ where the conductance 
through the dot followed the CST power law l|3"§jl . From 
the exponent, the interaction strength was deduced to be 
g = 0.54, corresponding to a TLL parameter g p = 0.23 
for the charged sector of the effective four-channel TLL 
theory of SWNTsA 9 - The quantum dot was formed by 
two nearby buckles acting as tunneling barriers. Since 
CST effects are predicted only for quite transparent bar- 
riers, it is instructive to estimate A. This is simpler for 
a single buckle used in earlier experimentsjii where also 
a power-law linear conductance G\b(T) was observed. 
The single-barrier case is analytically solvable^ and for 
a high barrier, 



Gib _ ^ 5/2 r(l 



^{A/DfinT/D) 2 ^- 2 . 



4G 2T(l/2 + l/ S ) 

The measurementsii for Gib yield together with g = 0.54 
a barrier transparency A ss 45 meV, taking a bandwidth 
of D = 0.5 eV; see Ref. 0- Assuming that the buckles 
have similar features when two are designed in a row, 
we can now establish a connection to the double-barrier 
case. In Ref. [lJa, e — 38 meV was measured, yielding 
A w 1.2e, consistent with our conclusions above. For 
CST to be operative, one needs finite level spacing, low 
(but not too low) temperatures, and not too high barri- 
ers. These conditions apparently were met in the SWNT 
experiments in Ref. I la Let us then comment on the 
Monte Carlo results of Ref. l28l where also the CST power 
law H39fl has been found. For a direct comparison, we de- 
termine A for the potential strength Ui mp used in Fig. 4 
of Ref. |2^. Using Eq. |JBJ, for the simulation parame- 
ters g = 0.6,C/ imp = 0.2D, and E s = ttD/20, we find 
A ~ 3.3e, again consistent with our conclusions. We note 
in passing that Eq. I|39|) has been obtained in Ref. 
starting from finite-range interactions among the elec- 
trons. The divergence has been regularized by summing 
up a selection of higher-order terms. However, this selec- 
tion was too strict j21 leading to the incorrect conclusion 
that finite-range interactions would be a prerequisite for 
CST to occur. As shown here, also zero-range interac- 
tions suffice, as long as the tunneling barriers are not too 
high. 

We emphasize that for Fermi liquid leads (g = 1) and 
T < e, one finds G max oc T _1 both within a UST 
and a CST analysis. Several researchers recently ap- 
proached the double-barrier TLL problem by considering 
weak interactions, g close to 1, without evidence for CST 
scalingiS^SSiSS As we have discussed in Sec. IIVI in this 
weak-interaction limit, the crossover temperature T c very 
quickly goes to zero when decreasing A, excluding CST 
effects for weak interactions even for finite level spac- 
ing. Put differently, for g close to 1, the master equa- 
tion approach for large A will always break down (£ be- 
comes larger than 1) upon lowering T before CST sets in. 
Therefore our results are in fact consistent with previous 
results Furthermore, the functional renormaliza- 
tion group calculation of Ref. 29 reported traces of an "ap- 



parent" power law (as opposed to true scaling) that could 
be linked to the CST mechanism. Finally, Komnik and 
GogolirJi presented an exact solution of a related model 
at the point g = 1/2. In their model, however, there is no 
sequential tunneling regime at all, and therefore we be- 
lieve that this represents a non generic situation that has 
nothing to say about the issues at stake here. This point 
has also been clarified in other recent publicationsi 2 ^ 2 ^ 
Unfortunately, this also excludes the possibility of an in- 
dependent analytical check of our results. 

To conclude, we hope that these novel features of in- 
teracting one-dimensional electrons will stimulate other 
theoretical work as well as further experimental checks 
of the CST versus UST picture put forward here. 
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APPENDIX A: FOURIER EXPANSION SCHEME 



To evaluate the rate expressions in Sec. IIIII we often 
need the quantities e ±Wdot ^ and e ±w& ^\ Exploiting the 
periodicity of these correlation functions with period r e = 
2ir/e, it is convenient to perform a Fourier expansion, 



,Wdot(t) 



-w dot (t) 



-w A (t) 



,w A (t) 



= Y c p( e ) e 

p— — OO 
OO 

= Y dp(e)e 



P— — OQ 



= Y v P^ e 

p— — oo 
oo 

= Y w p( £ ) e 



-ipet 



-ipet 



-ipet 



-ipet 



P— — OQ 



We note that Eqs. (|23|) and (|27f) imply that these Fourier 
coefficients are real. For low temperatures T, but keeping 
leading corrections in y = e _£ / T , one finds, with the 
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Heaviside function 6(x): Note that codo = vqWq — 1. The Fourier coefficients are 

shown in Figs. HUIandllllfor g = 0.6 and T«£. In that 



c k {e,T) = c k - -y (c fe _ x + c k+l + 2c k ) + 0{y' z ) 



case, the T = form of the correlation functions 1 -?•!!) and 



g (|27[) can be taken. Since the number of non zero Fourier 

coefficients is not exceedingly large, a quick, vt 
rate, and reliable numerical scheme can be imp 
for the evaluation of diagrams of the first class. 



, (1 — e~ £ / D )~ 1 / 9 coefficients is not exceedingly large, a quick, very accu- 

Cfe = 8(k)[ — 1) — rate, and reliable numerical scheme can be implemented 



r(i/g + i) kE/D 

m/g-k + l) 



1 



dk(s,T) — dk H — y (dk-i + d k +i — 2dk) + 0(y 2 ), Next, as mentioned in Sec. IIIIBl we show in detail 

9 how diagrams of the first class (ft = 1, 2, 7, 8) in TableHH 

dk = 6{k) — — 6 e ^ D ) 1 ' /g £(Vg + ^) c -ks/D have been handled, taking r^'^(n) as concrete example. 

ft! r(l/g) We consider the expressions on resonance and for low 

Vk(s T) = (— l) k \Ck(s T) temperature T <C s. After inserting the above Fourier 

_ e i D ^_il g expansions, the integration over 72 can be performed di- 

v _ L jl e j rectly. Since wo«o = 1, the '— 1' in the square bracket 

^ ' expression in 7 , see Table Ull is exactly canceled by 

w k (e,T) = (—l) k x~ 1 dk(£,T). the corresponding Fourier term. We find 



A 4 ^ 

16 ^-^ ■« — ' 7 

fe,i=0 m,r,p,q=0 



dud, 



k aiw m v r v p w q 



IE- 



mr—pq JQ 



/>oo />oo 

/ d 7ie - w °>>»< ri > +fe -*-"« ,ri / dT 3 e" 

JQ JQ 



(Al) 



where we have introduced the notation e_ mr _ P9 = (— to+ such that 
r — p + q)e, and analogously for E- k - mp . The prime in 

indicates that the term with m = r = p = q = 0is 
excluded from the sum. Next, we define the dimension- 
less kernels [cf. Eqs. and lf25|) ] 

/>oo 

K K {E) = e / dTe~ Soh '^ cos[Et - R ohm (T)} 
Jo 

= ^rohzn(S), (A2) 

/>OC 

#i(£) = £ / dre- So ^ sin [Br - i?ohm(r)] , 
Jo 



r./.(2)/ \ o ^ 4 ^ V*' d k diw m v r v p w q . 
r fl,V W = - 2 77T^ 2^ 2^ 2 , ,2 — {7^R(-£fcmp)^R(e-i- Pg ) -7^i(-efcmp)^i(e-z- Pg ) 

fc,2=0 m,r,p,«=0 7 £ -rm--pg 

& —rnr—pq [Kn(-ekmp)Ki(e-i-p q ) + Ki{—E kmp )Kft,{E-i-pq)]} ■ 

I 

Since one needs to self-consistently enforce 7 <C £, we can with the 7-independent terms 
simplify this expression to the form 

r £,7 2) W = ^ (~^ 9 ,7 + ^ + C fl) r) , (A3) 
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Ag,7 = gE E dkdid r -p+ g c r Cpdq{Kn(—ek r q)Kji(e-i-p q ) — Ki(—ekrq)Ki(e-i-pq)} 

fc,i=0 r,p,g=0 
oo 00 

B ^ = -\ E E" w^ca {gR( _ e )gR(e ) 

8 ^-^ ' (— m + r — p + o) 

k,l=0 m.r,p,q=0 y 1 

-Ki(-e km p)Ki(e-i-pq)} 
C *S = A E E" dkdl * mCrCpd * (-ir+*-+P+* {-K^-e^K^e+pq) 

k,l—0 m,r,p,q—0 

-Ki(-E kmp )K n (E-l-pq)} . (A4) 



Here, the second prime in denotes the sum with the 
constraint —m + r — p + q 0. The remaining Fourier 
sums are performed numerically. 

For completeness, we finally summarize the corre- 



sponding results for the remaining diagrams of this type, 
k = 1,2,8. The final result for each diagram will be of 
the form l|A3|) . with A g! i = A g ^ = and 



1 °° °° i 

A?3 = o E E dkdic r - p+q d r dpC q {K- R (-e k rq)K K {s-i p -q) + Kl(-e k rq)Kl(e-l p - q )} , 



B 9- 



8 

fc,Z=0 r,p,g=0 



oo 00 

B ^ = 4 E E' J^f^zh (-l) m+r+P+q {Knhe kmp )K R (e^pq) 



fc,Z=0 m,r,p,q=0 y y H > 

-K\{-e kmp )Ki{e-i pq )} , 



#3,2 - 



sE E' r^rS^ (-D"' + - +p+g {^(-^K(- e(H ) 



' fm + r + n + o) 

fc,Z=0 m,r,p,q=0 y y H ' 

-Ki(-e kmp )Ki(-ei pq )} , 



OO 00 

^E E" rZ^™^ {K R (-e kmp )K R (s^ q ) 



' f— m + r — p + q) 
+ifi(-e femp )i ; i:i(e_;p_ 9 )} , 

OO 00 

^ = 4E E' dkd f m " rCpdq (-1)™+^ {-K R (-e kmp ) Kl (e- lpq ) 

fc,Z=0 m,r,p,g=0 ^ y 

-lfi(-efe mp )i ; ir R (e_ ; pg)} , 

OO 00 

C 9 , 2 = -gE E' tTr m tn + l (-l) m+ " +P+g {^R(- efa np)gl(- e « M ) 
fc,i=0 ro,r,p,g=0 ^ 

+/fi(-efe mp )i ; ir R (-e / p 9 )} , 

-i^l(-efcmp)-^R(e-ip-g)} • (A5) 



Finally, we explain why we expect C g (T) to be T in- 1,2,7,8). Using Eqs. 1|A4|) and (|A5|) . the only T- 
dependent, focusing on diagrams of the first class (fc = dependent terms are the K R /i(E) defined above. Al- 
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though Kr(E = 0) has a power-law T dependence with 
K R (E = 0) cx T 1 /a- 1 i for finite energies E > T, and 
hence for T <C £, the T dependence is exponentially sup- 
pressed since the finite level spacing acts as an effective 
energy bias, see Eqs. IJA3I) and Similar arguments 



apply to Kr(E), which can be evaluated numerically in 
a straightforward manner. Although no closed analytical 
expression can be given, the overall T dependence follows 
directly from these considerations. 



1 S. Tomonaga, Prog. Theor. Phys. 5, 544 (1950). 

2 J.M. Luttinger, J. Math. Phys. 4, 1154 (1963). 

3 K. Schonhammer, J. Phys.: Condens. Matter 14, 12783 

(2002) . 

4 T. Giamarchi, Quantum Physics in One Dimension 
(Clarendon, Oxford 2004). 

5 O.M. Auslaender, A. Yacoby, R. de Picciotto, K. W. Bald- 
win, L. N. Pfeiffer, and K. W. West, Phys. Rev. Lett. 84, 
1764 (2000). 

6 A.M. Chang, Rev. Mod. Phys. 75, 1449 (2003). 

7 C. Dekker, Phys. Today 52 (5), 22 (1999). 

8 R. Egger and A.O. Gogolin, Phys. Rev. Lett. 79, 5082 
(1997); Eur. Phys. J. B 3, 281 (1998). 

9 C. Kane, L. Balents, and M.P.A. Fisher, Phys. Rev. Lett. 
79, 5086 (1997). 

10 M. Bockrath, D.H. Cobden, J. Lu, A. G. Rinzler, R. E. 
Smalley, L. Balents, and P. McEuen, Nature (London) 
397, 598 (1999). 

11 Z. Yao, H.W.Ch. Postma, L. Balents, and C. Dekker, Na- 
ture (London) 402, 273 (1999). 

12 J. Nygard, D.H. Cobden, M. Bockrath, P.L. McEuen, and 
PE. Lindelof, Appl. Phys. A 69, 297 (1999). 

13 H.W.Ch. Postma, M. de Jonge, Z. Yao, and C. Dekker, 
Phys. Rev. B 62, R10653 (2000). 

14 C.L. Kane and M.P.A. Fisher, Phys. Rev. B 46, 15233 
(1992). 

15 H.W.Ch. Postma, T. Teepen, Z. Yao, M. Grifoni, and C. 
Dekker, Science 293, 76 (2001). 

16 D. Bozovic, M. Bockrath, J.H. Hafner, C. M. Lieber, H. 
Park, and M. Tinkham, Appl. Phys. Lett. 78, 3693 (2001). 

17 I.O. Kulik and R.I. Shekhter, Sov. Phys. JETP 4, 308 
(1975). 

18 C.W.J. Beenakker, Phys. Rev. B 44, 1646 (1991). 

19 G.L. Ingold and Yu.V. Nazarov, in Single Charge Tunnel- 
ing, NATO ASI Ser. B 294, ed. by H. Grabert and M.H. 
Devoret (Plenum, New York 1992). 

20 A. Furusaki and N. Nagaosa, Phys. Rev. B 47, 3827 (1993). 

21 M. Sassetti, F. Napoli, and U. Weiss, Phys. Rev. B 52, 
11213 (1995). 

22 A. Furusaki, Phys. Rev. B 57, 7141 (1998). 

23 A. Braggio, M. Grifoni, M. Sassetti, and F. Napoli, Euro- 
phys. Lett. 50, 236 (2000). 

24 M. Thorwart, M. Grifoni, G. Cuniberti, H.W.Ch. Postma, 
and C. Dekker, Phys. Rev. Lett. 89, 196402 (2002). 

25 Yu.V. Nazarov and L.I. Glazman, Phys. Rev. Lett. 91, 
126804 (2003). 

26 D.G. Polyakov and I.V. Gornyi, Phys. Rev. B 68, 035421 

(2003) . 

27 A. Komnik and A.O. Gogolin, Phys. Rev. Lett. 90, 246403 
(2003). 

28 S. Hiigle and R. Egger, Europhys. Lett. 66, 565 (2004). 

29 V. Meden, T. Enss, S. Andergassen, W. Metzner, and K. 
Schonhammer, Phys. Rev. B 71, 041302 (R) (2005). 

30 U. Weiss, Quantum Dissipative Systems, 2nd edition 



31 
32 



33 
34 
35 



(World Scientific, Singapore, 1999). 

Y. Hu and S. Mukamel, J. Chem. Phys. 91, 6973 (1989). 
R. Kubo, M. Toda, and N. Hashitsume, Statistical Physics 
II, Springer Series in Solid State Physics, Vol. 31 (Springer, 
Berlin, 1985). 

V. Weisskopf and E. Wigner, Z. Phys. 63, 54 (1930). 
R. Zwanzig, Physica (Amsterdam) 30, 1109 (1964). 
W.H. Press et al, Numerical Recipes in FORTRAN (Cam- 
bridge University Press, Cambridge, England 1992). 
H. Baur, A. Fubini, and U. Weiss, Phys. Rev. B 70, 024302 
(2004). 
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summed up. 
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(N,N') 

c 

X • 

N-2, N-2 

b 

x • x , 

N-1.N-1 c 

a 

x • X 

N, N o 



x a' 

FIG. 1: Relevant part of the (N, jV')-plane of the RDM 
(schematic). Diagonal states are indicated by filled circles, off- 
diagonal states are marked by crosses. We use the shorthand 
notation a = (N, N+l), b = (N-l, N), and c = {N-2, N-l), 
and complex conjugate states are indicated by the prime. 
For the irreducible Fq contribution to the N — » N — 1 rate, 
we have four jumps. One starts from (N, N) and ends in 
(N — 1,N — 1), visiting an intermediate diagonal state after 
every second jump. 




N, N 



FIG. 2: (Color online) Paths in the RDM for (a) uncorrelated 
sequential tunneling (UST), and (b) cotunneling (COT). Wig- 
gled lines schematically indicate "bath-induced" correlations 
for first (red)- and second (blue)-order transitions in To. Dia- 
gram (a) involves two irreducible golden rule transition rates, 
i.e., there are no correlations across the intermediate diagonal 
state. Diagram (b) is not considered in what follows since we 
study a conductance peak. 
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(a) indirect 
k=7 




(a) 



N-2, N-2 



(b) direct: CST1 




^3 N-2, N-2 



>x 2 

^1 N-1,N-1 



N, N 



FIG. 3: (Color online) Examples for CST diagrams of order 
Tq oc A 4 : (a) the "indirect" diagram k = 7 in Tables Q and |TT| 
and (b) the "direct" diagram CST1. 
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FIG. 4: Temperature dependence of the dimensionless pa- 
rameter C g in Eq. I|29^ for various g and D = lOe. Dashed 
lines are guides to the eye only. Notice the double-logarithmic 
scales. 



FIG. 5: (Color online) (a) Linewidth 7 and (b) the parameter 
£ in Eq. (1.'->:it as a function of T for different g at A = 6e, D = 
lOe. Dashed lines in (a) represent 4r (1) cx T~ 1+1/fl . Notice 
the double-logarithmic scales. (Using these parameters, the 
master equation approach breaks down for g = 0.9.) 




FIG. 6: (Color online) Same as Fig.QJa), but for A = 3e. 





T/e 



FIG. 8: (Color online) Temperature dependence of the con- 
ductance maximum G max (solid lines) for A = 3e and TLL 
parameters g = 0.6 and g — 0.7. Note the double-logarithmic 
scales. Dotted (dashed) lines represent the CST (UST) power 
law G max oc T qc st (Gmax x r^T). l nS et: G max for 
A = O.le and g between g = 0.6 (bottom) and g = 1.0 (top). 
The slopes coincide with oust — —2 + 1/g. 
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FIG. 10: Fourier components d k and c k (inset) for g = 0.6 
and T < e. 
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FIG. 11: Same as Fig. HOI but for the Fourier components Wk 
and Vk (inset). 



